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Abstract 

A microscopic calculation of the perpendicular current in doped multiple 
quantum wells is presented. Interface roughness is shown to affect the reso- 
nant transitions as well as to cause a nonresonant background current. The 
theoretical characteristics exhibit several branches due to the formation of 
electric field domains in quantitative agreement with experimental data. 

I. INTRODUCTION 

The electric transport in semiconductor superlattices is usually dominated by resonances 
between the localized energy levels inside the wells resulting in peaks in the current-field re- 
lation. This may yield complicated current-voltage characteristics exhibiting may branches 
due to the formation of domains of different electric field inside the sample (see, e.g., |l] 
and references therein). These experimental features could be qualitatively reproduced by 
theoretical models combining rate equations for the transport between the wells and Pois- 
son's equation While these approaches assume interwell transition rates which are 

either fitted or obtained from phenomenological models, we have recently proposed a way 
to calculate the transport microscopically [|J. We obtained good quantitative agreement 
with the experimental data of Ref. |5j for highly doped samples, where the scattering from 
ionized impurities causes a strong broadening of the levels. Here we consider the lower doped 
sample used in consisting of a 40 period GaAs/AlAs superlattice (barrier width 6 = 4 
nm, well width w = 9 nm, period d — b + w, doping Np = 1.5 x 10 11 /cm 2 per well, cross 
section A = 1.13 x 10 -4 cm 2 ) and investigate the impact of interface roughness which both 
contributes to the broadening and causes nonresonant transitions between the wells. 

II. THE BASIC TRANSPORT MODEL 

In the case of weakly coupled quantum wells the appropriate basis set is a product 
of Wannier functions ty u (z — nd) of subband v localized in well n, and plane waves e*--. 
Here the z direction is the growth direction and k, r are vectors within the (x, y) plane. 
Restricting ourselves to the lowest two minibands (denoted by a and b) and coupling between 
neighbouring wells the Hamiltonian in the presence of an electric field F is given by H = 
H + Hi , where 
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Ho = J2 [K(k)aUk)a n (k) + E b n (k)bi(k)b n (k) - eF(J#4 (k)b n (k) + i#&t O^K®) 
Hi = £ [^ a < +1 (A;)a n (A;) + T*bt +1 (k)b n (k) - eFRfal +x {k)b n {k) - eF R\ a bl +1 {k)a n {k) 



with E^(k) = E u + h 2 k 2 /(2m w ) — eFdn [m w is the effective mass in the well), the couplings 
R v h v = I dz^ u (z — hd)z^/ u (z), and the miniband width 4|Tf| of subband za Diagonalizing 
the Hamiltonian Hq leads to renormalized coefficients in Hq and H\ Q which we use in the 
following. 

We calculate the Wannier functions in a Kronig-Penney-type model. Following Ref. || 
we model the nonparabolicity by an energy dependent effective mass m(E) = m c (l + (E — 
E c )/E g ), where m c is the effective mass at the conduction band minimum of energy E c , 
and Eg is the energy gap. Then the usual connection rules hold for the envelope function 
provided that the momentum matrix element P between the conduction and valence band 
states is identical in both materials. We use the values |8| = 0.067m e , m A1As = 0.15m e , 

^GaAs = L52 eVj £AiAs = 3.13 eV, and the conduction band discontinuity AE C = 1.06 eV. 
These parameters yield a relation E(k) = E c + h 2 k 2 / (2m(E)) which is in excellent agreement 
with the band structure of AlAs for the energies of interest .Q We obtain the coefficients 
E a = 47.1 meV, E b = 176.6 meV, 7\ a = -0.0201 meV, = 0.0776 meV, R\ a = 2.66 x 10" 4 rf, 
and R b a = -0.149d. 

For small couplings between the wells and fast intersubband relaxation the current from 



subband v in well n to subband fi in well n + 1 is given by the following expression Jl0| : 

r°° df 

I V nZ» + i = 2e E \H^ +mnM \ 2 j ^ ^KiM, S) ■ K + x(t, £ ~ A/i) [n F (S) - n F (S - A//)] . (1) 

Here A^(k, £) is the spectral function of subband v in well number n and np{£) = (l+e^) _1 
is the Fermi function. The energy £ is measured with respect to the electrochemical potential 
/i n in well n yielding an implicit dependence of A£(k,8) on /z n . We determine /i n from 
the local electron density n n = Y!,kS?5>od£TiF(£)A% l (k,£)/('7rA). Then the difference A// = 
/i n +i — A*n is equal to —eFd for n n = n n+ i. We obtain A^{k,£) in equilibrium from the 
retarded self-energy S^ ret (fc, £ ) neglecting the coupling to the other wells. In Ref. [§J we have 
calculated the self-energy for scattering from the screened potential of ionized impurities 
within the self-consistent single-site-approximation. As an additional contribution to the 
self-energy we study here the impact of interface roughness. 



III. MODELLING OF INTERFACE ROUGHNESS 

We consider an interface located at z = zq exhibiting thickness fluctuations £(r) of the 
order of ±i] (we use rj = 2.8 A which is one monolayer of GaAs). We assume (£(r)) r = 



1 Within this approach P = E g j '(2m c ) is slightly different for the two materials in contrast to 
the assumption. Furthermore, the envelope functions for different energies are not orthogonal as 
the effective Hamiltonian is energy dependent. However, the overlap is small and we neglect these 
complications. 
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and the correlation 



(C(l)f(r'))r = ar/ 2 exp(-|r-r'|/A). (2) 

This can be motivated in the following way: At a given point r there is an island of thickness r\ 
with a probability a (we use an average coverage a = 0.5). Therefore (£(r)£(r)) r = otrf. Pro- 
vided the island extends from r to r' we assume a constant probability to find a further neigh- 
bouring atom beyond r' yielding the exponential in Eq. (0). Following Ref. (Tl| we model 
the additional potential by a ^-function at the perfect interface U(r, z) = £,(r)AE c 5(z — z ) 
and obtain 



trough = }_ ^ [^(g)^^ + R )a n (k) + U» b (p)bl +h (k + p)b n {k) 

+U b h a (p)bl +h (k + p)a n (k) + U a h b (p)al +h (k + p)b n (k)\ (3) 



with the matrix elements U^(p) = f d 2 re-^ r -AE c [£(r)ty u *(z - hd)V(z )]. The elements 
contribute to the current from one well to the next via Eq. (p]). For weakly coupled 
wells U% M , U^, . . . are small and are neglected in the following. The elements result 
in scattering within the wells. We calculate their contribution to the self-energy within the 
self- consistent Born- Appr oximat ion 

yaret /. f] = 2 V I ^0° (& ~ jj) 1 2 m 

rough U, ^ ^2 £ + £F _ Ek , _ 5) ' 1 ^ 

where the factor 2 takes into account the two interfaces per well. The calculation for the 
subband b is performed in the same way. With ([|) we obtain the matrix element 

\U-{P)? = AAE^Mzo)\ 4 -^^ 2 (5) 

(1 + (pA) 2 ) 



In contrast to this the usual choice (see, e.g., [[Tl|]) (£(n)£(r')) r = arj 2 exp(—\r — r'\ 2 / X 2 ) 
yields \UQ a (p)\ 2 ~ exp(— A 2 p 2 /4). In order to compare both expressions we consider the 
matrix element for scattering at a single circular island of thickness rji and radius Aj which is 
given by U^ at (p) = 27ir]iAE c \^/ a (z )\ 2 \iJi(p\i) / p, where J\ is the Bessel function of the first 
kind. Now there are oiA/{ti\ 2 ) islands and we assume that negative and positive rji cancel 
each other. Then \U aa {p)\ 2 = Y,i \Uo at (p)\ 2 ~ P 3 f° r l ar g e V as m Eq. (M with a prefactor 
differing by 2/ it. This indicates that the choice (0) captures an essential part of physics not 
contained in the usual choice. 

Furthermore, the single island result reflects the quality of the Born-approximation for 
the self-energy. The analysis of a 5-potential reveals that diagrams containing multiple 
scattering at a single island become important if the product U(p)m/(2h 2 ) is larger than 
1. Using the value UQ ai (0) this yields 2h 2 /(irm\ 2 ) < r]AE c \^/ a (z )\ 2 . (For our parameters 
A > 15 nm.) This expression can be interpreted easily: The right hand side is the energy the 
electron gains if it is located at an island exhibiting a larger well width. The left hand side 
is the quantization energy associated with a length A. The localization of the wavefunction 
takes place if the energy gain due to the larger well width dominates the cost due to the 
localization. In this case higher order diagrams become important. 
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IV. RESULTS AND DISCUSSION 



Summing the contributions between different subbands given by Eq. (M we obtain the 
current 7„_ >n+ i(i 7 ', n n , n n+ i) depending on the field and the electron densities. Fig. la shows 
I n ^ n+ i(F, Np, Nd) exhibiting a first maximum I = 0.33 mA at eFd = 10 meV where the 
current is dominated by T" transitions and a second peak I = 1.8 mA at eFd = 140 meV 
due to the a — > b resonance. The interface roughness has two implications: Firstly the 
maxima are broadened and slightly shifted due to the contribution to the self-energies via 
Eq. (£|). Secondly there is a nonresonant current due to the Ui a transitions dominating the 
behaviour between the maxima. 

In order to study the domain formation we consider effective fields F{ between wells i 
and i + 1 fulfilling the discretized Poisson equation e r eo(Fi — i^_i) = e(rij — No)- Then the 
local currents I^i + i(Fi, m, have to be equal for all i in the steady state. The total 
voltage across the superlattice is U — dJ2iLi F%- Using the boundary conditions n\ = 2N D , 
n^+i = 2Nd we obtain the current- volt age characteristic shown in Fig. lb which exhibits 
the usual sequence of branches . In contrast to previous theoretical results the maximal 
current of the branches (90 //A) is significantly lower than the height of the first resonance. It 
is almost independent of the boundary conditions, which mainly affect the voltage where the 
leftmost homogeneous branch breaks up. The maximal current of 90 /iA is between the two 
experimental results 60 /iA and 130 /iA for samples #1 and #2, respectively, which both 
have the nominal sample parameters used here. We can reproduce the difference in current 
between both samples by assuming barrier widths which differ by just one monolayer, thus 
suggesting a possible explanation for the experimental discrepancy. The minimal currents of 
the branches (26 /iA) seem to be lower than the experimental values. They depend strongly 
on the nonresonant currents, e.g., we obtain the value 10 /iA ignoring transitions via U\. 
This indicates that we have probably underestimated the nonresonant current. The slope 
of the branches is steeper in our calculation than found experimentally; this could be due 
to some additional contact resistance of the order of 1 kQ which is not included in our 
calculation. 

In conclusion, we have presented a microscopic model for sequential tunnelling in weakly 
coupled quantum wells. The model describes the domain formation in superlattices quan- 
titatively with no fitting parameters. Nonresonant transitions due to interface roughness 
have been shown to strongly affect the current between the resonances and the extent of the 
current branches. 
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FIGURES 
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FIG. 1. (a) Current-field relation for equal densities in neighboured wells. The full and dotted 
lines are calculated with and without interface roughness, respectively (assuming A = 7 nm). (b) 
Current-voltage characteristics for a superlattice consisting of N = 40 wells for voltage sweep-up. 
The inset gives the field profile at U = 2.5 V. 
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